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ABSTRACT 

Using three-dimensional cosmological simulations, we study the assembly process of one of the 
first galaxies, with a total mass of ~ 10 8 Mq, collapsing at z ~ 10. Our main goal is to trace the 
transport of the heavy chemical elements produced and dispersed by a pair-instability supernova 
exploding in one of the minihalo progenitors. To this extent, we incorporate an efficient algorithm 
into our smoothed particle hydrodynamics code which approximately models turbulent mixing as a 
diffusion process. We study this mixing with and without the radiative feedback from Population III 
(Pop III) stars that subsequently form in neighboring minihalos. Our simulations allow us to 
constrain the initial conditions for second-generation star formation, within the first galaxy itself, and 
inside of minihalos that virialize after the supernova explosion. We find that most minihalos remain 
unscathed by ionizing radiation or the supernova remnant, while some are substantially photoheated 
and enriched to supercritical levels, likely resulting in the formation of low-mass Pop III or even 
Population II (Pop II) stars. At the center of the newly formed galaxy, ~ 10 5 Mq of cold, dense 
gas uniformly enriched to ~ 10~ 3 Zq are in a state of collapse, suggesting that a cluster of Pop II 
stars will form. The first galaxies, as may be detected by the James Webb Space Telescope, would 
therefore already contain stellar populations familiar from lower redshifts. 
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1. INTRODUCTION 

One of the most important goals in modern astro- 
physics is to understand the end of the cosmic dark ages, 
when the first stars and galaxies transformed the simple 
early universe into a sta t e of ever increasing complexity 
(IBromm fc Larsod 12001 iCiardi fc Ferraral |2005t iGloverl 
120051 iBarkana fc Loebl I2007t iBromm et al.l 120091 ). The 
first stars, the so-called Population III (Pop III), were the 
source of hydrogen-ionizing UV photons, thus initiating 
the e xtended process of cosmic reionization (Fa n et al.1 
2006) . They also synthesized the first heavy chemical el- 
ements, beyond the hydrogen and helium produced in the 
big bang, to be dispersed into the pristine intergalactic 
medium (I GM) through supernova (SN) explosion s and 
winds (e.g.. iMadau et al.ll2001t iWise fc Abelll2008D . 

An intriguing possibility for the first stars is that some 
of them may have died as a pair-instability supernova 
(PISN), a peculiar fate predicted for prog enitor masses 
in th e range 140 Mp, < M* < 260 M (jBarkat et al.l 
119671 iHeger fc Wooslevl 120021 ) . Current theory pro- 
poses a top-heavy initial mass function (IMF) for the 
first stars, with a characteristic mass > 100 Mq 
(lAbel et alJl2002t IBromm et alJl200lfO'Shea fc Normanl 
120071: lYoshida et al.ll2008D . Together with the expecta- 
tion that mass loss due to ra diatively-driven winds is 
negligible at low metallicities (Kudritzki 2002|), one ar- 
rives at the robust expectation that at least a fraction of 
Pop III stars should have died as PISNe. Compared to 

tgreif@mpa-garching.mpg.de 



conventional core-collapse SNe, a PISN is distinguished 
by not leaving behind a compact remnant. Instead, the 
exploding star is completely disrupted, and all metals 
produced inside the Pop III star are released into the sur- 
roundings, leading to me tal yields of y = M z /M* ~ 0.5 
(jHeger fc Wooslevi [2001 ) . An abundant occurrence of 
PISNe in the early universe could thus have rapidly es- 
tablished a bedr ock of meta l s, at least locally, of order 
Z > 10~ 3 Z Q (|Greif et al.l [20071) . Recently, the ex- 
tremely luminous SN 2007b i was observed in a nearby 
galaxy (|Gal-Yam et al.l [20090 . with characteristics, such 
as very large Ni masses, that seem to unambiguously 
point to a PISN origin. This discovery greatly strength- 
ens the possibility for finding these events at high red- 
shifts as well. 

We here carry out cosmological simulations tracing the 
detailed assembly process of a primordial galaxy, tak- 
ing into account the feedback effects from Pop III star 
formation inside the mini halo progeni t or sy stems. Our 
work extends the study bv lGreif et al.l (|2008l ). which fol- 
lowed the virialization of gas in the galactic potential 
well under the idealized assumption of no such feed- 
back. Specifically, we include radiative feedback, leading 
to the build-up of H II regions around Pop III stars (e.g., 
Whalen et all 120041: lAlvarez et al.l 120061 : lYoshida et ail 
2007aHGreif et al.l l2009b) , as well as the mechanical and 



chemical feedback from a single PISN that explodes in 
the earliest minihalo progenitor. The latter feedback 
refers to the additional cooling that becomes available in 
metal-enriched gas, allowing the gas to reach lower tern- 
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peratures, and to possibly fragment into low-mass Pop- 
ulation II (Pop II) stars. It has been suggested that this 
transition in the star formation mode, from high-mass 
Pop III to normal-IMF Pop II, occurred once a mini- 
mum metal enrichment was in place, th e so-called criti- 
cal metallicity Z r . rit ~ 10~ 6 -10~ 35 Z Q (jBromm fc Loebl 
l2003t I Schneider et al .1 120 06) . It is therefore important to 
identify any "super-critical" regions within a simulation 
box to distinguish between Pop III and Pop II star for- 
mation sites. The overall goal is to predict the proper- 
ties of the hrst galaxies, to be observed with the James 
Webb Space Telescope fJW ST), planned for launch in 
2014 ([Gardner et all [20061. One key ingredient neces- 
sary for such predictions is to quantify the amount and 
distribution of metals inside the emerging galaxy, just 
prior to the onset of its initial starburst. This would 
allow us to constrain the detailed mix of stellar popula- 
tions, and consequently to arrive at observables such as 
broad-band colors, emission line signatures, and lumi- 
nosities (e.g.. Uohnson et aT1l2009|) . 

A complementary empirical probe of the first stars 
is given by "Stellar Archaeology", the study of abun- 
dance patterns in metal-poor Galactic halo stars 
(|Beers fc Christliebl [2001 : IFrebel et al.ll2007L l200l . and 
of glo bular clusters in the Milky Way and other galaxies 
(e.g., iHelmi et al.ll2006t iPuzia et al.f[2006h . To infer the 
characteristics of the first SNe, and therefore to constrain 
the primordial IMF from the observed fossil chemical 
record, we need to better understand the in situ physical 
conditions in the regions where the first low-mass stars 
formed. This again requires us to simulate the transport 
of metals from the first SNe in a realistic cosmological 
setting. A related problem is the apparent absence of a 
PISN abundance patte rn in any of the known metal-poor 
stars ([Tumlinsonl 12006) . It has been argued that such a 
PISN signature may be hidden in stars with relatively 
high metallicities, Z ~ 10~ 2 5 Zq, due t o the extremely 
high yield from even a single explosion (Karlsson et al. 
2008). Our simulations can assess the viability of this 
scenario with much improved physical realism. 

We trace the mixing and dispersal of metals, origi- 
nating in a PISN inside a Pop III minihalo at z ~ 30, 
all the way to their reassembly into the growing poten- 
tial well of the first galaxy at z ~ 10. The physics 
governing this transport is highly complex. Some as- 
pects, such as the advection of metals during hierar- 
chical structure formation, is reliably modeled by our 
Lagrangian smoothed-particle hydrodynamics (SPH) ap- 
proach. Other processes, such as hydrodynamical in- 
stabilities and turbulence, are not well resolved in our 
simulation. We may therefore miss the fragmentation 
of gas within the dense shell of the SN remnant, which 
could result in the formation of an intermediate gen- 
eration of stars (IMackev et alj|2003t IQ'Shea fc Normanl 
l2lMIWhaTen et al.ll2008bl ; iNagakura et al.ll2009D . How- 
ever, we have developed an efficient algorithm to approxi- 
mately model the mixing that r esults from these su b-grid 
effects as a diffusion process (IGreif et all l2009ah . The 
corresponding diffusion coefficient is evaluated assum- 
ing a simple mixing-length approa ch to turbulent trans - 
port even in the supersonic regime ()Klessen fc Linl[200l . 
with quantities that are entirely local. This allows us 
to greatly improve on our earlier, "ballistic " transport 
of discrete metal packets (|Greif et al.l 120071) . and con- 



sequently to derive meaningful metallicity distribution 
functions. 

The structure of our work is as follows. In Section 2, 
we describe our numerical methodology and the proce- 
dure to set up the initial conditions. We then discuss 
the assembly process of the first galaxy, and the result- 
ing distribution of metals within our simulation volume 
(Section 3). In Section 4, we investigate the gas prop- 
erties in the different sites where second-generation star 
formation could occur. We conclude by summarizing and 
assessing the cosmological implications of our results. All 
distances quoted in this paper are in proper units, unless 
noted otherwise. 

2. NUMERICAL METHODOLOGY 

2.1. Simulation setup 

We perform our simulations in a cosmological box of 
size 1 Mpc 3 (comoving), initialized at z — 99 with a fluc- 
tuation power spectrum corresponding to a concordance 
A cold dark matter (ACDM) cosmology with matter den- 
sity ft m = 1 — 51a = 0.3, baryon density fit, = 0.04, Hub- 
ble parameter h = f/o/100 km s _1 Mpc^ 1 = 0.7, where 
Hq is the present Hubble expansion rate, spectral index 
n s = 1.0, and a normalization erg = 0.9 ([Spergel et al.l 
2003) . We take the chemical evolution of the gas into ac- 
count by following the abundances of H, H + , H , H2 , H^~ , 
He, He + , He ++ , and e~, as well as the three deuterium 
species D, D + , and HD. We include all relevant cooling 
mechanisms, i.e., H and He collisional ionisation, excita- 
tion and recombination cooling, bremsstrahlung, inverse 
Compton co oling, and collisional exc itation cooling via 
H 2 and HD ([Glover fc Jappsenl l2007V We explicitly in- 
clude H2 cooling via collisions with protons and elec- 
trons, which is important for the chemical and thermal 
evolution of gas with a significant fractional ionization 
(|Glover fc Abel|[200l . 

In a preliminary run with 64 3 gas and dark matter 
(DM) particles, we locate the formation site of the first 
~ 10 8 Mq halo. This object is massive enough to cool 
through Lya emission and fulfi l our prescription for a 
"first galaxy" ([Greif et a l. 2008). We subsequently carry 
out a standard hierarchical zoom-in procedure to achieve 
high mass resolution in the comoving volume of the 
galaxy (e.g..lNavarro fc Whitd H994: Tor men et al1li"997l: 
IGao et al.ll2005D . We apply three consecutive levels of re- 
finement, such that a single parent particle is replaced by 
a maximum of 512 child particles. The particle masses 
for DM and gas in the highest resolution region contain- 
ing the comoving volume of the galaxy with an extent 
of ~ 300 kpc (comoving) are ~ 33 M and ~ 5 M Q , 
respectively. The mass resolution is ~ 400 M , so that 
the Jeans mass at a density of njj = 100 cm~ 3 and tem- 
peratures near the cooling threshold set by the cosmic 
microwave background (CMB) is well resolved. 

2.2. Feedback by the first star 

The first minihalo collapses at a redshift z ~ 30 near 
the location of the nascent galaxy. Once its central den- 
sity exceeds nn = 100 cm~ 3 , we assume that a sin- 
gle Pop III star with 200 Mq forms, near the mid- 
dle of plausible values de termined in numerical studies 
(|Q'Shea fc Normanl [20071) . However, for the purpose of 
determining the properties of the ensuing H 11 region, 
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we use the ionization, diss ociation, and h eating rates 
for a 50 M Q Pop III star (jSchaererl 120021) . which lim- 
its the severity of the radiative feedback. We note that 
the SN remnant will likely propagate somewhat further 
in this case, since pressure equilibrium with the ambi- 
ent medium will be achieved at later times (see also 
iGreif et "ail 120071 ). On the other hand, the photoheat- 
ing from nearby stars may mitigate this effect and lead 
to only minor deviations from the true metal distribu- 
tion. The corresponding surface temp erature, luminos - 
ity, and stellar lifetime may be found in lSchaererl (2002), 
while detail s of the ray-tracing algorithm are described in 
IGreif etldl {2009b) ■ Two key assumptions of this model 
are that the gas is optically thin to molecule-dissociating 
radiation, and that the effects of stellar evolution are 
negligible, which would otherwise lead to a decrease in 
the number of i onizing photons emitted at the end of the 
main sequence ([Marigo et al.|[2001t ISchaereri r2002). 

At the end of its brief lifetime, the progenitor is pre- 
dicted to explode as a PISN, ejecting o f order 100 M B in 
meta ls and 10 52 erg of kinetic energy (|Heger fe Wooslevl 
2002). In the simulation, we assign the above energy 
and metal content to t he cen tral 200 Mq of gas as de- 
scribed in IGreif et al.l (|2007| ). The explosion is initial- 
ized at the end of the free expansion phase, where the 
momentum of the swept-up mass leads to a transition 
to the Sedov-Taylor phase. This simplified treatment is 
roughly valid as long as the central density is low, which 
is est ablished by the previ ous photoheating from the star 
(e.g.. I Johnson et al.ll2007h . All particles within the ini- 
tial ejecta are assigned a metal yield of 50 %, consisting 
of C, O, and Si and relative abundances corresponding 
to a solar composition. Alth ough this is somewhat in - 
consistent with the results of iHeger fc Wooslevl (2002), 
who show that the actual C and Si abundances differ 
roughly by a factor of four from the solar values, the 
oxygen abundance is reproduced correctly. Since this is 
the most common element, we expect this caveat to not 
qualitatively affect our conclusions. We note that other 
elements such as Mg, S, Ca, and Fe are produced as well, 
but are negligible in ter ms of providing additional cooling 
(ISantoro fc ShulII [2006h . 

2.3. Mixing of metals 

A key difficulty in modeling chemical mixing in SPH 
simulations is that no inherent mass flux between parti- 
cles exists. Most simulations simply work around this 
cavea t (e.g., IScannapieco et all 120051 : IKobavashi et al.l 
12007ft , while others implement ad hoc methods to smoot h 
the metals within the kernel (e.g.. iWiersma et al.ll2009l ). 
A more realistic treatment clearly requires a physical 
model for the distribution of metals between particles. 
We have recently devel oped a diffusion-ba sed algorithm 
to accomplish this task (Greif ct al. 2009a). On all scales 
of interest, chemical mixing proceeds via turbulent mo- 
tions that cascade down to smaller and smaller scales. By 
quantifying the degree of turbulence on the smallest re- 
solved scale, one can estimate the local degree of mixing. 
The exchange of metals between particles is then mod- 
eled as a diffusion process, with the diffusion coefficient 
set by the smoothing length and the velocity dispersion 
within the kern el. The SPH equ ations for diffusion are 
well known (see Monaghan 2005), and only minor modi- 
fications are required to apply this machinery to chemi- 



cal mixing. We note that a similar method has recently 
been use d in the context o f IGM enrichment by massive 
galaxies (|Shen et al.ll2009D . 

2.4. Metal cooling 

As a consequence of a more accurate metallicity esti- 
mate, we can model the additional cooling provided by 
heavy elements. For this purpose we assume that C, O, 
and Si are produced with solar relative abundances. This 
is a relatively good approximation for the cooling pro- 
vided by the heavy elements ejected by a 200 M Q PISN. 
For these species, our model distinguishes between two 
temperature regimes. At temperatures below 2 x 10 4 K, 
we use a slightly modified version o f the c hemical net- 
work presented in lGlover fc Jappsen (2007). This treat- 
ment follows the chemistry of C, C + , O, + , Si, Si + , 
and Si ++ , in addition to the primordial species men- 
tioned above, and accurately models the effects of fine 
structure cooling from C, C + , O, Si, and Si + . It does 
not include the effects of cooling from molecules such 
as CO or H2O, whic h is unimportant at the densities 
probed in this s tudy dGlover fc Jappsenll2007fi . A mod- 
ification to the iGlover fc Jappsenl (|2007| ) treatment is 
the use of more accurate rate coefficients for the colli- 
sional excitation of the fine-structure lines of C and O 
by atomic hydrogen, taken f rom the recent calculations 
of IAbr ahamss on"et al.l (|2007|) . The second, more impor- 
tant modification is that we keep C and Si in their first 
ionized states, since they are easily ionized by a soft UV 
background, which is likely to be present at z < 20 (e.g., 
IBromm fc Loeb|[200l . 

At temperatures above 2 x 10 4 K, a full non-equilibrium 
treatment of metal chemistry, such as we use in cold 
gas, becomes significantly more costly, owing to the in- 
creasing number of ionization states that become avail- 
able. We theref ore adopt a cooling rat e deriv ed from the 
values given in ISutherland fc Dopital (Tl99l for gas m 
collisi onal ionization equilibrium. ISutherland fc D opita 
(1993) give total cooling rates (including hydrogen and 
helium line cooling, and bremsstrahlung) for metallic- 
ities Z = 0, 10~ 3 , 10" 20 , lO" 15 , 10" 1 - , 10~ - 5 , 1.0, 
and 10 5 Zq. By subtracting off the values for the zero 
metallicity case from the other cases, we can isolate the 
contribution of the metals. For metallicities in the range 
10~ 3 < Z < 10° 5 Zq, we can then determine the metal 
cooling rate by interpolation of these values. For metal- 
licities Z < 10~ 3 Zq and Z > 10° 1,5 Zq, we must extrap- 
olate from these values, which is inevitably more inac- 
curate, although we note that in the former case metal 
cooling is ineffective compared to H and He cooling. 

In view of the significant uncertainties that ex- 
ist regarding the nature and the quantity of the 
dust produced by Pop III SNe (ISalvaterra et all 
2004t iBianchi fc Schneider! I2007t iNozawa et al.l 120071: 
Cherchneff fc D wck 2009|), we do not include dust in our 
present simulation. We do not expect this to lead to 
a significant error: although dust cooling may be ex- 
tremely important during the late stages of protostellar 
collapse in low metalli city gas (see, e.g., lOmukai et al.l 
120051 : IClark et aLl|2008 f). it is not an effective coolant at 
the gas densities resolved in this study. 

2.5. Feedback by second-generation stars 
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The thermal and chemical evolution of the IGM is com- 
plicated by the fact that minihalos in the vicinity of the 
SN remnant collapse and form stars, which photoheat 
the gas that is eventually incorporated into the galaxy. 
Additional feedback by accretion from stellar remnant 
black holes or subsequent SN explosions may become 
important as well. However, in the present work, we 
choose to not include more than one SN, since we be- 
lieve it is essential to understand the effects of a single 
SN before attempting to follow the evolution of multi- 
ple, interacting SNe. We therefore assume that all stars 
formed in minihalos collapse directly to black holes after 
they die. Recent investigations have shown that accre- 
tion onto these stellar remnant black holes is inefficient as 
a result of the low d ensity of the gas after pho toheating 
(j Johnson fc Bromml [20071 1 Alvarez et al.ll2009f ). at least 
for regions of the universe t hat are not part of an unusu- 
ally large overdensity (e.g-. IGao et al.ll2005|) . This leads 
us to neglect radiative feedb ack from early miniquasars 
(e.g., IKuhlen fc Madaul[200l . 

In light of the complexity induced by multiple star for- 
mation sites, we have chosen to run two distinct sim- 
ulations: in simulation A (Sim A), we insert a 50 Mq 
Pop III star whenever the density in a minihalo exceeds 
7iH = 100 cm -3 , and subsequently solve for the H II re- 
gion according to the prescription in Section 2.2. After 
of order 50 stars have formed, we stop the ray-tracing 
routine and instead form sink particles. This ensures 
that the simulation does not become exceedingly com- 
plex. Once a sink particle is created, all gas particles 
within the smoothing length corresponding to the above 
density threshold are immediately accreted. Further ac- 
cretion is governed by the criterion that a particle must 
fall within this smoothing le ngth (details on the i mple- 
mentation may be found in Uappsen et all ((2005) ) . In 
simulation B (Sim B), we always form sink particles and 
do not include any radiative feedback beyond that of the 
very first star. This distinction allows us to investigate 
the effects of photoheating on the distribution of metals 
and the assembly of the galaxy. 

3. PRE-GALACTIC ENRICHMENT 

In this section, we discuss the evolution of the SN rem- 
nant and its influence on the surrounding medium, fol- 
lowed by a detailed investigation of the distribution of 
metals in both simulations. This will then allow us to 
discuss the consequences of metal enrichment for second- 
generation star formation in Section 4. 

3.1. Assembly of the first galaxy 

The expansion of the SN remnant into the IGM is de- 
picted in Figures 1 and 2, where we show the central 
100 kpc (comoving) of Sim A and B. After ~ 15 Myr, 
the SN remnant has propagated to ~ 1 kpc in radius, 
roughly the size of the original H II region. Both simu- 
lations are identical at this point, since no other miniha- 
los have collapsed yet. In a few cases, the shock ahead 
of the SN remnant has disrupted neighboring miniha- 
los that are just in the process of collapsing. However, 
most minihalos remain unaffected since they are cither 
too far away or have collap sed to high enough densi- 
ties (|Cen fc Riquclmc 2008). The distinction between 
disrupted and largely unaffected minihalos leads to an 
interesting result that will be discussed in Section 4. 



After roughly 100 Myr, or a Hubble time at z ~ 30, 
the SN remnant comes into pressure equilibrium with 
the surrounding medium and stalls. During this pe- 
riod, multiple neighboring minihalos undergo runaway 
collapse and form stars. In Sim A, the ensuing H n re- 
gions then partially overlap with the SN remnant and the 
comoving volume of the nascent galaxy. At this point, 
the complexity of the simulation dramatically increases 
and a multi-phase medium forms. In Figure 3, we show 
the distribution of the gas in density and temperature 
space within the high resolution region at z ~ 10, ex- 
cluding gas accreted onto sink particles. It occupies more 
than six orders of magnitude in density, ranging from 
riH = 4 x 10~ 5 to 100 cm -3 , and more than three or- 
ders of magnitude in temperature, ranging from T = 10 
to 4 x 10 4 K. Most of the gas is shock-heated to the 
virial temperature and cools once the density exceeds 
7iH ~ 1 cm -3 , but some parcels of gas in Sim A are pho- 
toheated to 10 4 K and have cooled down to T ~ 1000 K 
by z ~ 10. This warm, diffuse medium resides in the low- 
density IGM at uh ~ 10~ 3 - 10~ 4 cm~ 3 and will not be 
available for further star formation for > 10 8 yr, the ap- 
proximate IGM free-fall time. A third phase, consisting 
of gas shock-heated by the SN remnant, has cooled down 
from an initial temperature of ~ 10 s K to ~ 10 4 K and 
has become indistinguishable from the relic H n region 
gas. 

Concomitant to the onset of star formation in other 
minihalos, the potential well of the nascent galaxy as- 
sembles at the center of a massive filament (see bottom 
panels of Figures 1 and 2). Accreted halos are subject 
to tidal stripping and rapidly lose their mass, similar to 
the disruption of satellite galaxies entering the halo of 
the Milky Way. In some cases, they complete a few or- 
bits before being disrupted, or undergo a major merger. 
Interestingly, virial shocks around halos in Sim B are 
hotter and more distinct than in Sim A. The additional 
photoheating in Sim A increases the pressure at larger 
radii, reducing the infall velocity and, in consequence, 
the post-shock temperature of the gas. For the same rea- 
son, structures in Sim B are generally more pronounced, 
which is particularly true at the center of the galaxy (see 
inlays). 

3.2. Distribution of metals 

During the first few million years after the SN ex- 
plosion, the metals are distributed into the IGM by 
the bulk motion of the blast wave, which leads to a 
tight correlation between temperature and metallicity. 
This becomes evident from the top panels shown in Fig- 
ures 1 and 2. Once the SN remnant stalls, mixing is 
facilitated by agents that act on various scales. On 
the largest scales, filamentary accreti on dominates (e.g., 
IWise fc Abeil \mm IGreif et all l2008h . while on smaller 
scales photoheating by neighboring stars as well as dy- 
namical interactions associated with the virialization of 
the galaxy become important. As can be seen in the mid- 
dle and bottom panels of Figures 1 and 2, these all inter- 
act to create a complex hierarchy of turbulent motions 
that efficiently mix the gas. In Figure 4, we show the 
resulting density-metallicity relation for the three out- 
put times in Figures 1 and 2. At early times, a distinct 
correlation arises: Hot, underdense regions of the IGM 
tend to be enriched to Z ~ 10~ 3 Z or higher, while the 
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Fig. 1. — From left to right: the density-squared weighted average of the hydrogen density, temperature, and metallicity along the line 
of sight in the central 100 kpc (comoving) of Sim A. From top to bottom: a time series showing the simulation 15, 100, and 300 Myr after 
the SN explosion. The inlays show the central 10 kpc (comoving) of the nascent galaxy. The metals are initially distributed by the bulk 
motion of the SN remnant, and later by turbulent motions induced by photoheating from other stars and the virialization of the galaxy. 
As can be seen in the inlays, the gas within the newly formed galaxy is highly enriched. 
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Sim B Width: 100 kpc (comoving) From top to bottom: 

Z in „ = 30 Inlays: 10 kpc (comoving) At = 15, 100, 300 Myr 




J 
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Fig. 2. — Same as in Figure 1, but without additional photoheating from nearby stars (Sim B). This results in a significantly smaller 
enriched region as well as less efficient mixing. The gas within the galaxy is more concentrated, and we find that virialization shocks are 
hotter and more pronounced than in Sim A, where the gas shocks at larger distances from the halo, so that the infall velocity is smaller. 
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Fig. 3. — Mass- weighted distribution of gas in density and tem- 
perature space within the high resolution region at z ~ 10, color- 
coded according to its mass fraction. Most of the gas is shock- 
heated to the virial temperature and cools once the density exceeds 
~ 1 cm -3 , but a significant fraction of the gas in Sim A is also relic 
H II region gas that has cooled to T ~ 1000 K and resides in the 
low-density IGM at n H ~ 10~ 3 - 10~ 4 cm" 3 . 

densest regions, such as nearby minihalos, remain almost 
devoid of metals. Once the potential well of the galaxy 
assembles, metal-rich gas accumulates at high densities 
and leads to a flattening of the relation. 

A more detailed view of the metal distribution is shown 
in Figures 5 and 6, where we plot the enriched mass 
as a function of gas overdensity and metallicity. Early 
on, the initial ejecta have not propagated very far and 
the metals occupy a small region in density and metal- 
licity space, located at the relic H II region density of 
~ 0.1 cm~ 3 and at supersolar metallicities. Over time, 
the SN remnant then distributes its metals to the sur- 
rounding medium and the average metallicity decreases, 
while the total enriched mass increases. Once the poten- 
tial well of the galaxy has become deep enough, a large 
fraction of the hot, metal-rich gas residing in the IGM 
has recollapsed to high densities. By z ~ 10, the en- 
riched mass contained within the densest parcels of gas 
in the galaxy has increased to roughly 10 5 Mq, with a 
clear peak in the distribution at Z ~ 10~ 3 Z@. Interest- 
ingly, this exceeds any critical metallicity quoted in the 
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Fig. 4. — Mass-weighted, average metallicity as a function of gas 
density within the high resolution region, shown ~ 15 (solid line), 
100 (dotted line), and 300 Myr (dashed line) after the SN explosion 
(Sim A). At early times, a distinct correlation is evident: under- 
dense region are highly enriched, while overdense region remain 
largely pristine. Once the potential well of the galaxy assembles, 
metal-rich gas becomes dense and the relation flattens (see also 
Figures 5 and 6). 

literature (|Bromm fc Loebl[200l ISchneider et "all [20061. 
and we expect that a stellar cluster with a normal IMF 
will form (e.g.JOmukai et alj|2005t ISchneider et alj|2006t 
IClark et al.ll2008l 12009ft . It therefore appears that a sin- 
gle PISN is sufficient to trigger a transition from Pop III 
to Pop II star formation. 

Finally, we note that the additional photoheating in 
Sim A has an interesting effect on the distribution of 
metals. As is evident from Figures 1 and 2, the en- 
riched region is generally larger and better mixed in 
Sim A. In some cases, the photoheating even ejects sig- 
nificant amounts of enriched gas out of the potential 
well of the galaxy, which can be seen in the case of the 
metal bubble extending to the lower left of the simula- 
tion box. This behavior is intriguing in light of claims 
that the IGM is substan t ially enriched at hig h redshifts 
(jSongaila fc Cowid 119961: ISchave et all 120031) . although 
this is generally attributed to SN feedback instead of 
photoheat ing, which acts on objects with small circula r 
velocities (|Thoul fc Weinberdll996l:lDiikstra et~aTll2004D . 
Comparing Figures 5 and 6, we also find more quantita- 
tive evidence of this effect. From top left to bottom right, 
an increasing amount of enriched gas accumulates at low 
densities, which then accretes onto the nascent galaxy. 

3.3. Metal Cooling 

In both simulations, we have included fine-structure 
cooling for C, O and Si at low temperatures and colli- 
sional excitation and recombination cooling at high tem- 
peratures. Depending on metallicity, the latter can dom- 
inate over H and He cooling at very early times. Such a 
phase exists during the first few million years after the 
SN explosion, when the temperature of the SN remnant 
exceeds 10 5 K and the metallicity within the remnant 
is well above solar. By performing test runs with and 
without metal line cooling in a smaller simulation box, 
we have found that metals temporarily enhance the net 
cooling rate by a factor of a few. However, even for the 
case of a PISN, the influence on the dynamical evolution 
of the SN remnant remains small. 

At temperatures below 10 4 K, fine-structure cooling 
provided by heavy elements takes over. However, we 
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Sim A: Metallicity distribution function 
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Fig. 5. — Enriched mass as a function of metallicity and gas overdensity, shown for different times after the SN explosion (Sim A). The 
overdensity is color-coded according to the legend shown in the top left panel. At early times, the metals remain confined to the initial 
stellar ejecta, located at a density of ~ 0.1 cm -3 and at supersolar metallicity. Later on, the gas mixes with the surrounding medium 
and the total enriched mass increases, while the average metallicity decreases. Once the potential well of the galaxy assembles, the gas 
rccollapses to high densities and becomes available for star formation. At the final output time, the metallicity distribution within the 
galaxy peaks at Z ~ 10~ 3 Zq. 
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A Sim = Sim A — Sim B: Metallicity distribution function 
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Fig. 6. — Same as in Figure 5, except that we show the differential enriched mass obtained by subtracting the metal content in each 
metallicity bin in Sim B from that in Sim A. From top left to bottom right, the amount of enriched gas ejected into the IGM by the 
additional photoheating in Sim A increases. At the final output time, a large fraction of this gas has accreted onto the galaxy. 
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find that once the gas has cooled to the regime where 
fine-structure cooling becomes important, the metallic- 
ity has dropped such that cooling is instead dom i nated 
by H 2 . This was already found by Uappsen et al.1 (|2007l 
l2009a| ). who argued that dust is responsible for chang- 
ing the cooling and fragmentation properties of the gas. 
Since cooling via dust grains kicks in at densities much 
higher than we can resolve here, a detailed treatment is 
beyond the scope of this work, although we expect that 
it will pl ay a crucial role for the further evolution of the 
gas (e.g.. lOmukai et aT1l2005l ). We note that the relative 
importance of the different cooling channels could be al- 
tered in the presence of a strong Lyman- Werner (LW) 
radiation backg round, which may l ead to the dissocia- 
tion of H 2 (e.g. lHaiman et al.|[l997t Uohnson et~aT1l2007l; 
lAhn et al.ll2009h . 

4. SECOND-GENERATION STAR FORMATION 

To elucidate the pathway to the first galaxies, an es- 
sential ingredient is to characterize the conditions in the 
gas clouds out of which the s econd generation of stars 
form (e.g.. lBromm et aLll2009l ). We in particular need to 
quantify the metal content in such clouds. This applies 
to the first galaxy itself, but also to neighboring mini- 
halos which collapse during the assembly process. It is 
generally believed that minihalos can only host Pop III 
stars, as self-enrichment is unlikely, and transport of met- 
als out of the I GM into pre-condensed halos might be 
severely limited (|Cen fc Riquelmd 120081) . Here, we can 
address these questions within their proper cosmologi- 
cal setting. We note that we use the terminology to 
classify early stellar populations, in particular the dis- 
tin ction between Pop III .l and III. 2, recently introduced 
bv iMcKee fc Tan! (|2008t ). The former refers to metal- 
free stars whose formation is entirely unaffected by any 
stellar feedback, and is instead governed by cosmological 
initial conditions, whereas the latter denotes metal-free 
stars that form under the influence of feedback effects. 

4.1. Star formation in minihalos 

Concomitant to the assembly of the galaxy, nearby 
minihalos undergo runaway collapse and form stars. De- 
pending on distance and central density, their collapse 
could be unaffected, delayed, or completely inhibited 
in the presence of an ionization front or the SN rem- 
nant. Recent detail ed radiation-h y drodyn a mics si mula- 
tions performed by iWhalen etail (|2008al |R 12009ft have 
revealed that for a large region in parameter space 
most minihalos remain unaffected. However, in dis- 
tinct cases the ionization front or the SN shock can 
disrupt the core of the halo, heating and ionizing the 
gas and possibly leading to metal enrichment. This in 
turn may trigge r the formation of less massive so-called 
Pop I II. 2 stars (lO'Shea et al.l[2005t Uohnson fc Bromml 
120061: lYoshida et all l2007bl ). or normal Pop II stars if 
enough metals are brought to the center of the halo. The 
amount of fragmentation sensi tively depends on the ini- 
tial conditions of the collapse (jJappsen et al.ll2"009bt l. 

How does this compare with our results? In Figure 7, 
we show all star-forming minihalos in the high resolution 
region as a function of redshift and distance from the 
center of the galaxy. Crosses denote minihalos that re- 
main largely unaffected by feedback effects and will likely 
form Pop III.l stars. Minihalos that have been severely 
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Fig. 7. — All star-forming minihalos in the high resolution re- 
gion as a function of redshift and distance from the nascent galaxy 
(Sim A). Crosses denote undisturbed halos that will form Pop III.l 
stars, star symbols denote halos that have been photoionized and 
will likely form Pop III. 2 stars, and triangles denote halos that 
have been enriched to Z > 10 Zq by the SN remnant and will 
likely host normal Pop II stars. Most minihalos have collapsed to 
high enough densities by the time star formation ensues and thus 
remain pristine. 

disturbed by ionizing radiation from a nearby star show 
elevated electron and H 2 fractions and are denoted by 
star symbols. These halos also form significant amounts 
of HD and will most likely cool to the temperature of 
the CMB before becoming Jeans-unstable, following the 
canonical pathway of Pop III. 2 stars. Finally, triangles 
denote minihalos that have been disrupted by the SN 
remnant and are enriched to Z > 10~ 6 Z@. Two of these 
halos are even enriched to Z > 10 -3 ' 5 Z Q , and will al- 
most certainly form Pop II stars. 

In Figure 8, we show a particularly interesting case 
that elucidates the importance of timing in regulating 
the strength of feedback. Here, a minihalo has been par- 
tially photo-heated by a nearby star to 10 4 K, and sub- 
sequently enriched with metals. In the process of its col- 
lapse, the core has fragmented into two distinct objects: 
one in which the molecule fraction has been significantly 
enhanced due to the photoionization, allowing the gas 
to cool to ~ 100 K, and another in which the cooling 
is unaltered from the standard minihalo pathway. In 
both cases, the metallicity is not high enough for metal 
fine-structure cooling to become important, although at 
higher densities dust-induced cooling may take over and 
lead to the formation of a small cluster of Pop II stars. 
This case suggests that the formation of low-mass stars 
need not await the formation of second-generation halos 
with M vir > 10 s M Q . 

4.2. Star formation in the nascent galaxy 

What are the properties of the gas within the newly 
virialized galaxy? From Figures 5 and 6, it appears that 
~ 10 5 Mq of cold, dense gas at its center have been en- 
riched to Z ~ 10 -3 Zq. This is confirmed by Figure 9, 
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Fig. 8. — Radial distribution of gas in a minihalo that has been 
partially shock-heated by an ionization front and subsequently dis- 
rupted by the SN remnant (Sim A), color-coded according to its 
mass fraction. The temperature panel shows that the core has 
fragmented into two distinct clumps that are both enriched to 
Z ~ 3 X 10~ 4 Zq: one in which the photoionization has led to 
enhanced molecule abundances, with H2 and HD fractions > 10 — 3 
and > 10 — 8 , respectively, allowing the gas to cool to 100 K, and 
another where the gas has shielded itself from radiation and has 
only cooled to ~ 200 K. The CMB temperature at these redshifts 
is ~ 50 K. Metal fine-structure cooling does not become impor- 
tant, although dust-induced cooling may set in at higher densities 
than we can resolve here. 




0.01 0.1 1 10 100 

n H [cm" 3 ] 



og M/M 



-5.0 -4.5 -4.0 -3.5 -3.0 -2.5 

Fig. 9. — Distribution of gas in density and temperature space 
within the virial radius of the galaxy (Sim A), color-coded accord- 
ing to its mass fraction. The dense gas at the center of the galaxy 
is enriched to Z ~ 10 — 3 Zq and will likely form of a stellar clus- 
ter with a more normal IMF. Similar to the minihalo case shown 
in Figure 8, the photoheated gas has formed enhanced molecule 
abundances and will likely cool to the temperature of the CMB at 
somewhat higher densities. 
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where we show its distribution in density and temper- 
ature space within the virial radius of the galaxy. It 
has become highly enriched by accretion from the sur- 
rounding IGM and is in a state of collapse. We also find 
an interesting effect related to photoheating: previously 
ionized gas forms elevated H2 and HD abundances and 
coexists with unheated gas. Similar to the minihalo case 
discussed in Figure 8, the gas might therefore later frag- 
ment and form individual clumps. We note that in the 
density range resolved here, cooling is dominated by pri- 
mordial molecules in s tead of metal fine-structure cooling 
(jJappsen et all 120071 l2009aT) . At somewhat higher den- 
sities, we expect that the gas will cool to the temper- 
ature of the CMB before the equation of state hardens 
again. However, a more detailed investigation of the sub- 
sequent fragmentation must await high-resolution simu- 
lations that resimulate the central few kpc of the galaxy. 

5. SUMMARY AND CONCLUSIONS 

We have performed a set of highly resolved SPH simu- 
lations that allow us to investigate the enrichment of the 
IGM by a PISN exploding in a high-redshift minihalo. 
We have incorporated a substantially higher degree of 
realism compared to our previous work in the form of 
radiation feedback, explicit chemical mixing, and metal 
line cooling. In particular, we have employed a physi- 
cally motivated model for the mixing of met als between 
individ ual SPH particles based on diffusion (|Greif et al.l 
I2009aj) . We have followed the distribution of metals as 
they are ejected into the IGM and then recollapse into the 
larger, M v i r ~ 10 8 M Q potential well of a "first galaxy" 
assembling at z ~ 10. We have performed simulations 
with and without radiative feedback to assess the effects 
of photoheating on the assembly of the galaxy and the 
distribution of metals. 

The heavy elements ejected by the SN are initially dis- 
tributed into the IGM by the bulk motion of the SN 
remnant, until it stalls and mixing is instead facilitated 
by turbulent motions on smaller scales. These are in- 
duced by the dynamics of the underlying DM, and, to a 
smaller degree, by the additional photoheating. A clear 
correlation between metallicity and gas overdensity is es- 
tablished, where the densest regions consisting of exist- 
ing minihalos remain largely pristine, while voids around 
the SN progenitor become highly enriched. This corre- 
lation breaks down once the potential well of the galaxy 
assembles and metal-rich gas residing in the IGM recol- 
lapses to high densities. The metallicity at the center of 
the galaxy then grows to Z ~ 10 -3 Zq, likely resulting in 
the formation of a stellar cluster with a more normal IMF 
(|Clark et al.ll2008fl . Although the influence of photoheat- 
ing on the mixing of the gas is limited, an interesting side 
effect is that it can expel enriched gas out of the potential 
well of the galaxy into the IGM. Finally, the mechanical 
feedback exerted by the photoheating and the propaga- 
tion of the SN remnant is usually quite small, although in 
some cases nearby minihalos are disrupted and enriched, 
possibly leading to the formation of Pop III. 2 and Pop II 
stars. 

One of the most promising observational tools to un- 
derstand the nature of these objects relies on the ele- 
mental composition of second-generation stars. Previ- 
ous searches for PISN signatures have targeted the most 



metal-poor stars in our Galaxy , with an iron abundance 
as low as [Fe/Efl ~ 10" 5 (e.g., Beers fc Christliebl [20051: 
iFrebel et al.l 120051) . Unfortunately, the spectra of these 
stars do not reveal the strong odd-even effect indicative 
of PISNc, implying that they were either not very com- 
mon, or that existing surveys did not target a suitable 
sample of stars. Indeed, our simulations show that a 
single PISN is sufficient to enrich a gas cloud quite uni- 
formly to Z ~ 10 -3 Zq, implying that its signature may 
be buried in stars with a significantly higher metallic- 
ity (jKarlsson et al.ll2008l) . The ideal search strategy may 
therefore have to be modified to concentrate on stars with 
[Fe/H] ~ 10 -3 . Recent surveys of extragalactic globular 
clusters lend support to this idea, since the oldest of these 
have been found to contain large [ q/Fe] ratios, which is 
indicative of enrichment by PISNe (jPuzia eTaUHool). 

The standard picture that minihalos only form Pop III 
stars, or remain sterile before being disrupted d uring sub- 
sequent merging fe.g.. lHaiman fc Hol der 2003), may also 
have to be modified. Whereas that is indeed the case 
for the majority of minihalos in our simulation, there 
exists a small number of halos that do contain dense, 
metal-enriched gas. Although we do not follow its sub- 
sequent evolution, it is intriguing to speculate whether 
these clouds might fragment to form the first globu- 
lar cl usters, or at least sligh t ly less massive pr ecursors 
(e.g.. IBromm fc Clarkd [20021: IBolev et~aL1 l2009h . Soon, 
it should become possible to study the further fate of 
these second-generation star forming clouds with the re- 
quired extremely high resolution by resimulating the cen- 
tral halo gas, and by including dust cooling and opacity 
effects at high density. Similar simulations will be able to 
address the formation of the first star clusters inside the 
first galaxies, starting from the initial conditions deter- 
mined here. The goal of understanding galaxy formation 
from first principles, at least in the relatively simple case 
of dwarf-sized systems at the beginning of hierarchical 
structure formation, might finally have come into reach. 
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